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Abstract 

We study effects of eccentricity fluctuations on the elliptic flow coefficient V2 at mid-rapidity in 
both Au+ Au and Cu+Cu collisions at y%5v = 200 Ge V by using a hybrid model that combines 
ideal hydrodynamics for space-time evolution of the quark gluon plasma phase and a hadronic 
transport model for the hadronic matter. We find that the effect of eccentricity fluctuation is 
modest in semicentral Au+Au collisions but significantly enhances V2 in Cu+Cu collisions. 



1. Introduction 

One of the major discoveries at Relativistic Heavy Ion Collider (RHIC) in Brookhaven Na- 
tional Laboratory is that the elliptic flow [Q]] is found to be as large as an ideal hydrodynamic 
prediction for the first time in relativistic heavy ion collisions |2|]. Systematic studies revealed, 
however, that the reasonable agreement between results from ideal hydrodynamics and elliptic 
flow data has been achieved only by a particular combination of dynamical modeling J3J]. If one 
looks at the comparison carefully, one finds the agreement between hydrodynamic results and the 
data is not perfect. This would be due to an absence of initial fluctuation effects 0]. Therefore, 
further investigation is indispensable toward better understanding of the elliptic flow data and the 
quark gluon plasma (QGP). 

2. Initial Conditions 

For initial conditions in this study, we employ the Monte-Carlo version of both the Glauber 
model [5] and the factorized Kharzeev-Levin-Nardi (fKLN) model [jf! to generate the initial 
distribution of entropy density on an event-by-event basis. We first calculate a transverse entropy 
density profile so(x_i_) = s(r = tq,x, y, r] s = 0) in each sample at an impact parameter b for a given 
centrality, where tq = 0.6 fm/c is the initial time for the hydrodynamical simulations. Then the 
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Here (• • ■) describes the average over transverse plane by weighting the entropy density in a 
single sample: 
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The eccentricity with respect to the reaction plane, the participant eccentricity, and the corre- 
sponding transverse areas can be defined 01, respectively, as 
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5rp = n^Jcrlcrj, S part = nJa^af^o%,. (6) 

The impact parameter vector and the true reaction plane are not known experimentally. So one 
can set an apparent frame of created matter shifted by (x,y) = ((x), (y)) and tilted by T* from the 
true frame in the transverse plane ]7[] : 

2 2 

tan2 v F= — -. (7) 

2<T xy 

The anisotropy of particle distribution could be correlated with the misidentified frame. To ac- 
count for this, we first shift the center-of-mass of the system to the origin in the calculation frame 
and then rotate the profile in the azimuthal direction by *P to match the apparent reaction plane to 
the true reaction plane. We generate the next sample of an entropy density profile again as above 
and average the profiles over many samples. We repeat the above procedure for many samples 
until the initial distribution is smooth enough. The initial conditions obtained in this way con- 
tain the effects of eccentricity fluctuation even though the profile is smooth. In particular, even in 
case of vanishing impact parameter, the eccentricity is finite due to its event-by-event fluctuation. 
It is the particle distribution calculated from the initial conditions mentioned above that can be 
directly compared with the experimental data. Note that the procedure of averaging over many 
samples without shift or rotation corresponds to conventional initialization without the effect of 
eccentricity fluctuation. 

By using the initial conditions above, we simulate space-time dynamics of matter produced 
in relativistic heavy ion collisions with a hybrid approach in which the macroscopic description 
of the QGP is followed by the microscopic description of the hadron gas. For details, see Ref. ■ 




Figure 1: Eccentricity as a function of iVpart in (a) Au+Au and (b) Cu+Cu collisions. Solid, dotted, dashed and dash- 
dotted lines correspond to e part from the MC-KLN model, £ palt from the MC-Glauber model, £rp from the MC-KLN 
model and £rp from the MC-Glauber model, respectively. 
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3. Results 



Figure Q] shows the eccentricities as functions of the number of participants N p!at with or 
without eccentricity fluctuations. The results from the MC -Glauber and MC-KLN models are 
compared with each other in (a) Au+Au and (b) Cu+Cu collisions. In semicentral Au+Au col- 
lisions at 10-50% centrality, the effect of initial fluctuations enhances the eccentricity parameter 
by 8- 1 1 % (5-8%) in the MC-Glauber (MC-KLN) model. The enhancement factor s pdn /s RP is the 
largest at the very central bin (0-5% centrality): £ pa rt/£RP = 1-83 in the MC-Glauber model and 
1.53 in the MC-KLN model. A qualitatively similar behavior is observed in Cu+Cu collisions, 
but quantitatively the effect of fluctuation is significant due to the smallness of the system. 
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Figure 2: Centrality dependences of V2 for charged particles at mid-rapidity in (a) Au+Au and (b) Cu+Cu collisions 
at ysAW = 200 GeV with the MC-Glauber model initial conditions are compared with PHOBOS data (filled plots) |9]. 
Open circles (squares) are results with (without) eccentricity fluctuation. 
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Figure 3: Centrality dependences of v% for charged particles at mid-rapidity in (a) Au+Au and (b) Cu+Cu collisions at 
■fitm = 200 GeV with the MC-KLN model initial conditions are compared with PHOBOS data (filled plots) @], Open 
circles (squares) are results with (without) eccentricity fluctuation. 

Figures [2] and [3] show the centrality dependence of V2 for charged hadrons at mid-rapidity 
(| J] |< 1) in Au+Au and Cu+Cu collisions at ysAW = 200 GeV. From Figs.|2](a) and (b), one 
sees that the effect of eccentricity fluctuations is not significant in Au+Au collisions, whereas 
V2 is largely enhanced in Cu+Cu collisions due to fluctuation effects. These tendencies are also 
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expected from the results of initial eccentricity as shown in Fig.Q] It is interesting to point out 
that elliptic flow coefficients \>2 in the Glauber model initial conditions still slightly undershoot 
the experimental data, particularly in Cu+Cu collisions, even with fluctuation effects. Thus there 
is almost no room for the viscosity in the QGP stage to play a role in the Glauber initial conditions 
within the hybrid approach with the ideal gas equation of state in the QGP phase. Figure[3]is the 
same as Fig. [2] but the initial conditions are taken from the MC-KLN model. Again, a similar 
qualitative behavior is seen: The effect of eccentricity fluctuations is small in Au+Au collisions 
but is large in Cu+Cu collisions. Due to larger initial eccentricity in the MC-KLN model than 
the MC-Glauber model, the results are somewhat larger than the experimental data in peripheral 
Au+Au collisions, whereas we reasonably reproduce the v>2 data with initial fluctuation effects 
in Cu+Cu collisions. Viscous effects in the QGP phase could reduce the v'2 and enable us to 
reproduce the data in Au+Au collisions in this case. However, the results are already comparable 
with the data in Cu+Cu collisions even though the number of participants is almost the same as 
that in peripheral Au+Au collisions. So it would be non-trivial to establish whether the same 
viscous effects also give the right amount of V2 in Cu+Cu collisions. 

4. Conclusions 

We calculated the elliptic flow coefficient as a function of the number of participants in the 
QGP hydro plus the hadronic cascade model and found that the effect of eccentricity fluctuations 
is visible in very central and peripheral Au+Au collisions and is quite large in Cu+Cu collisions. 
This strongly suggests that the effect of eccentricity fluctuations is an important factor which has 
to be included in the dynamical model for understanding of the elliptic flow data and for precise 
extraction of transport properties of the produced matter from the data. 

So far, one has been focusing on comparison of hydrodynamic results with V2 data only 
in Au+Au collisions. The experimental data in Cu+Cu collisions also have a strong power to 
constrain the dynamical models. Therefore, simultaneous analysis of V2 data in both Au+Au and 
Cu+Cu collisions will be called for in future hydrodynamic studies. 
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